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Abstract 

We use tools from g-calculus to study LDL T decomposition of the Vandermonde matrix V q with 
coefficients Vij = q lJ . We prove that the matrix L is given as a product of diagonal matrices and 
the lower triangular Toeplitz matrix T q with coefficients Uj = 1/0; q)i-j, where (z',q)k is the q- 
Pochhammer symbol. We investigate some properties of the matrix T q , in particular, we compute 
explicitly the inverse of this matrix. 
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1 Introduction and main results 


Let us consider a Vandermonde matrix 


K:= 


1111 
1 q q 2 q 3 
1 q 2 q A q 6 
1 q 3 q e q 9 


of size n x n. In the special case when q = e~ 2m l n , this matrix is called the discrete Fourier transform 
matrix. Explicit matrix factorizations of the discrete Fourier transform matrix are very important, since 
they are often used in various versions of the Fast Fourier Transform algorithm [3]. Motivated by this 
connection, in this note we plan to study the LDL 1 factorization of the matrix V q and to investigate the 
properties of the factors appearing in such decomposition. The tools and techniques, which are used to 
prove our results, come from g-calculus. 

First, let us present several definitions and notation. In what follows, we assume that n G N and 
q G C. We define the q -Pochharrirrier symbol 


A; q)n ■■= (1 - z){l - zq) ■ ■ ■ (1 - zq n 1 ), n> 1, 


( 1 ) 


*Dept. of Mathematics and Statistics, York University, 4700 Keelc Street, Toronto, ON, M3J 1P3, Canada. 
E-mail: kuznetsov@mathstat.yorku.ca 

Research supported by the Natural Sciences and Engineering Research Council of Canada. 


1 






and (z;q)o 1. We will denote by / the n x n identity matrix. The following matrices of size n x n 
will be used frequently in this paper: a lower-triangular Toeplitz matrix T q = {£jj}o<ij<n-i defined 
by tij = 1 /(q\q)i-j if i > j, and a diagonal matrix P q having coefficients {(< 7 ; q)i}o<i<n-i on the main 
diagonal, or, more explicitly, 


■ 1 

0 

0 

0 ..." 


"l 

0 

0 

0 

1 

1 

0 

0 ... 


0 

(t q ) i 

0 

0 

1 

(9:9)2 

1 

(9:9)1 

1 

0 ... 

P q := 

0 

0 

(<?; 9)2 

0 

1 

(9:9)3 

1 

(9:9)2 

1 

(9:9)1 

1 ... 


0 

0 

0 

(959)3 ••• 











Note that the matrices T q and T q - 1 are well-defined for all q G C \ A n , where the set A n is given by 

A n {q £ C : = 1 for some j — 1 , 2 ,..., n — 1 }. 

In our first result we identify explicitly the matrices appearing in the LDL T factorization of the 
Vandermonde matrix V q . 

Theorem 1. Assume that q e C \ A n . Then V q = LDL T , where L = P q T q (P q )- 1 and D is a diagonal 
matrix having coefficients {(—l)*g*b _ 1 h 2 (g; g)j} 0 <i< n -i on the main diagonal. 

In section 2 we give a very simple proof of Theorem 1 (our proof is based on the (/-Binomial Theorem). 
Alternatively, one could derive this result starting from formulas (2.4) and (2.5) in the paper [ 2 ] by Oruc 
and Phillips, who use symmetric functions to study LU decomposition of general Vandermonde matrices. 

In our second result we present some properties of the Toeplitz matrix T q , including an explicit 
formula for its inverse. First we define the following two matrices of size n x n: 


S := 


0 0 0 0 
10 0 0 
0 10 0 
0 0 10 


D q := 


10 0 0 
0 q 0 0 

0 0 q 2 0 
0 0 0 q 3 


( 2 ) 


Theorem 2. Assume that q e C \ A n . 

(1) (Tf - 1 = T q -i(I - S) — D q -iT q -iD q . 


(ii) Form EN we have 


T q D q -mT q -lD q n 


771—1 


'+E 

3 = 1 


(g 1-m ; q)i oi 

( 9 ; v)j 


(3) 


Remark 1. Note that the matrix H (I — S') -1 , which appears in item (i), is a lower triangular 
Toeplitz matrix having coefficients h t j — 1 if i > j and h,„j = 0 otherwise. Similarly, the matrix in the 
right-hand side of (3) is a lower-triangular Toeplitz matrix, having m non-zero diagonals: this matrix 
has coefficient 1 on the main diagonal and the coefficient (g 1-m ; q)j/(q] q)j on the sub-diagonal number 
j, for 1 < j < m — 1 . 
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2 Proofs 


The only tool that will be needed for proving Theorems 1 and 2 is the q-Binomial Theorem (see 
[1] [Theorem 10.2.1]), which states that 


jaz] q) op 

0 ;<?) oo 


E 


(g;g)/ ’ 


kl < i, kl < i. 


(4) 


Here (z\ q )^ := ~ z Q l ) and it is clear that this infinite product converges for all z G C and |g| < 1. 

We also record here the following two corollaries of the q-Binomial Theorem, which will be needed later: 


1 

(z]Q)oo 

0 ;g)oo 


£ 

j> o 


(<?; q). 


-, \q\ < 1 , \z\ < 1 , 


E -TT 3 -kl < 1 , zee. 


j> 0 


(?; q)i 


(5) 

( 6 ) 


Proof of Theorem 1: First we check that the matrix L = P q T q (P q ) 1 is a lower-triangular matrix with 
coefficients 

1 . = _hill_ i > , 

' J (« 9 ) 3 ( 9 ; 9 ) 4 - 3 ’ 

Considering an element (i, j) of the matrix LDL T we see that formula V q = LDL T is equivalent to the 
following identity: for any integers i,j >0 


min(ij) 

"= £ 


q lJ = 


T ) k q k{ - k 1)/2 (q-,q)i(q;q)j 


k=0 (p,q)k(q-,q)i-k(q-,q)j-k 

We will prove the above identity by writing the Taylor series of the function 

(■ uv ; q)oc 


g(u,v) := 


(u; q)oo(v, q) 

in two different ways. First of all, from formula (4) we obtain 


, \u | < 1 , h < 1 , kl < 1 , 


g(u,v) = 


x 


(uv, q )c 


(v, q )00 («; g)oo (v; 9)00 “ (g; g)* 


v 


(7) 


Using the fact that (v,q)i/(v;q) 00 
(5) we conclude that 


l/(g*u; g)oo and expanding this expression in Taylor series in v via 


g(u,v) = EE 

i >0 j >0 


q l 3y l v 3 

(g; g)*(g; g)j ’ 


( 8 ) 


On the other hand, we can obtain the series expansion of g(u,v ) by applying formulas (5) and ( 6 ) in 


3 















the form 


(uv;q) oo 

1 

(u;q) oo 

1 

{v,q) oo 


> y - -7-^— -W V 


k> 0 


E 

Z>0 

E 


m> 0 


(q\q)h 


u 


( 9 ; 9 ) 1 ’ 


( 9 ; 9)r 


We multiply the above three series expansions and obtain a Taylor series representation in the form 

_l^kgk(k- 1)/2 


9 K 77 ) = 


k> 0 Z>0 m>0 


(q;q)k(q;q )z(9; 9)» 


V+V +m 


(9) 


Comparing the coefficients in front of the term v l v 3 in both formulas ( 8 ) and (9) gives us the desired 
result (7). □ 


Proof of Theorem 2: Let us prove the identity T q T q -1 = (/ — S') -1 , which is equivalent to the first 
equality in item (i) (the second equality in (i) follows from formula (3) with m — 1). The main idea of 
the proof is that the Toeplitz matrix T q can be expressed in the following form 


T 

-L n 


/+ E 

3 > 1 


s j 

( 9 ; 9)/’ 


( 10 ) 


where the matrix S is defined in (2). The above formula is easy to derive, given that for 1 < j < n — 1 
the coefficients of the matrix S '- 7 have value 1 on the sub-diagonal number j and value zero everywhere 
else. In particular, S '- 7 is a zero matrix for j > n, thus the series in (10) terminates at j — n — 1. Similarly, 
using the identity 


( 1 / 9 ; 1 / 9 ); = (-1 ) J 9 j{j+ 1 )/ 2 (q;q)j. 


( 11 ) 


and formula ( 10 ) we obtain 


T q -i = I 


E 

3 > 1 


— 1 ) tgt 0'- 1 )/ 2 


( 9 ; q)i 


■(qsy 


( 12 ) 


Now, assume that |g| < 1. Then formulas (5) and (10) give us 


T q = [(S'; 9 )oo ] _1 = (/ - S)- 1 x (/ - qS)- 1 x (J - ^S )- 1 x • • • . (13) 


Similarly, formulas ( 6 ) and (12) give us 

T q - 1 = (qS- q)^ = (/ - qS) X (/ - q 2 S ) X (/ - g 3 S) X • • • . (14) 

From the above two identities we see that all the terms (/ — q l S) in the product T q T q -i are cancelled, 

except for the first term (/ — S') -1 , thus we obtain T q T q - 1 = (/ — S') -1 for |g| < 1. We extend this result 

from \q\ < 1 to the general case q E C \ A n by analytical continuation in q. 
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The proof of formula (3) uses the same ideas. Again, first we assume that |g| < 1. From (11) we 
check that D q -mT q -iD q m is a Toeplitz matrix of the form 


D q — mT q —l D q m — I ^ ^ 

3> 1 


JfiJti -!)/2 


-i ) j q- 


(?; q ) 3 




Using the above result and formula (13) we obtain 

T q D q -mT q -iD q m = [(^g)oo]- 1 x (q l ~ m S] q , ) 00 = (g^S; g) m _r. 
The desired desired result (3) follows by applying (4) and analytical continuation in q. 


□ 
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